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Abstract 



The classical and the quantal problem of a particle interacting in one- 
dimension with an external time-dependent quadratic potential and a con- 
stant inverse square potential is studied from the Lie-algebraic point of view. 
The integrability of this system is established by evaluating the exact invariant 
closely related to the Lewis and Riesenfeld invariant for the time-dependent 
harmonic oscillator. We study extensively the special and interesting case of 
a kicked quadratic potential from which we derive a new integrable, nonlin- 
ear, area preserving, two-dimensional map which may, for instance, be used 
in numerical algorithms that integrate the Calogero-Sutherland-Moser Hamil- 
tonian. The dynamics, both classical and quantal, is studied via the time- 
evolution operator which we evaluate using a recent method of integrating the 
quantum Liouville-Bloch equations The results show the exact one-to-one 
correspondence between the classical and the quantal dynamics. Our analysis 
also sheds light on the connection between properties of the SU(1, 1) algebra 
and that of simple dynamical systems. 
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I. INTRODUCTION 



The classical and quantum mechanical study of time- dependent Hamiltonian systems 
are generic and important. They span a wide spectrum of subjects ranging from interaction 
between atoms and radiation |2||| , adiabatic HH and non-adiabatic |||7| Berry phase, time- 
dependent harmonic oscillators [p|-[T0[1 and quantum motion of a particle in a Paul trap 
Ll], 12], to time-dependent mean- field theory (cranking model [13]) in particular and the 
time-dependent shell model [113 in general. 

The question of the existence of invariants (constants of motion) is one of central impor- 
tance in the study of any dynamical system, be it classical or quantal. A basic theorem of 
classical mechanics asserts that if the number of independent invariants, satisfying certain 
conditions, is equal to the number of degrees of freedom, then the motion can be reduced 
to quadratures, or equivalently, an action-angle transformation to a Hamiltonian dependent 



only on the actions can be found fl5|| . Such systems are integrable. 

The lack of sufficient number of invariants invariably leads to a phase space that has a 
nonzero measure of chaotic trajectories f!6 ]. For time-independent Hamiltonian systems the 
Hamiltonian itself is an invariant. However, when the Hamiltonian is an explicit function of 
time, it is no more an invariant, and this is of course a reflection of the non- conservation of 
energy. Various methods have been used to obtain approximate solutions for time-dependent 
problems, e.g., the adiabatic approximation, the sudden approximation, time-dependent 
perturbation techniques, etc. 

The most widely studied time-dependent Hamiltonian system is the time-dependent 
harmonic oscillator (TDHO). It has long been a problem of considerable interest because 
of its varied applications in different areas of physics, for instance in molecular physics, 
quantum chemistry, quantum optics and plasma physics. The Hamiltonian of this system is 
given by, 



n{t) = \v 2 + 



(i) 



where p and q are conjugate canonical variables. The adiabatic invariant for this system was 
originally given at the first Solvay Congress in 1911 when the Hamiltonian of this system 
was used as an approximate Hamiltonian for the slowly lengthening pendulum [17] . 

The study of this one- dimensional TDHO was greatly advanced due to the work of Lewis 

Lewis H determined the exact invariant by applying Kruskal's 



and Riesenfeld 10 



asymptotic method and showed that a previously known adiabatic invariant was in fact 
an exact invariant. Later Lewis and Riesenfeld (LR) [T(J determined that same invariant by 
starting with the assumption of the existence of an explicitly time-dependent, homogeneous 
and quadratic invariant of the form given by, 



I(t) = -[a(t)p 2 + (3(t)q 2 + 2 1 (t)pq}, 



(2) 



where the coefficients a(t),(3(t) and j(t) are time-dependent real functions and I(t) satisfies 
the condition, 



dl_ 
~dt 



m 

dt 



+ {I,H(t)} = 0. 



(3) 
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Here {*, *} denotes the usual Poisson bracket. From the above two equations and after some 
calculations they derived the exact invariant as, 



m = \ 



2 2 / -2 1 \ 2 o 

P V + [P + g - 2pppg 



(4) 



with p(t) satisfying the subsidiary condition 



p + ^(t)p-l = 0. (5) 
p d 

This more complicated, nonlinear, differential equation represented an advance due to the 
reason that any particular solution p(t) of the above equation would give the exact invariant 
I(t) for all initial conditions of p and q. 

The above technique actually predates the work of Lewis and Riesenfeld in a slightly 
different context, and we briefly point this out. We note that the equation of motion of 
the TDHO is, d 2 q/dt 2 + w 2 (t) q = 0. This is of the same form as the time- independent, 
one-dimensional Schrodinger equation, if we assume that t represents the spatial coordinate 
and q, the wave- function. According to an early work by Milne (1930) fl9| , it is possible to 
solve this Schrodinger equation using the knowledge of any particular solution of the Eq. 
after proper identifications. Therefore, any particular solution of the subsidiary condition 
can also give the exact solution for the TDHO. 

Lewis || attempted to give an interpretation of I(t) as the most general homogeneous 
quadratic invariant possible for the Hamiltonian of the one-dimensional TDHO. A more 
natural and physical interpretation has been suggested by Eliezer and Gray |2(J in terms 
of two-dimensional auxiliary motion, i.e., in terms of a two-dimensional uncoupled TDHO. 
They showed that the above subsidiary condition Eq. (|5p is the radial equation of motion for 
this two-dimensional system and the invariant I(t) is proportional to the conserved angular 
momentum of this auxiliary motion. Giinther and Leach |2lJ interpreted I(t) in terms of 



canonical transformations and under their transformation the invariant I(t) became the 
Hamiltonian of the one-dimensional time-independent harmonic oscillator of unit frequency. 

Besides these previous interpretations of I(t), we can interpret the form of I(t) chosen 
by Lewis and Riesenfeld (LR) [[3I|] from the Lie-algebraic point of view. The Hamiltonian of 
the one-dimensional TDHO is formed by the dynamical variables and |g 2 . These two 
dynamical variables together with pq are generators of the closed SU (1,1) algebra under 
the Poisson bracket operation and I(t) was chosen in as the linear combination of these 



generators. Generalizing the LR-invariant from the Lie-algebraic point of view, we expect 
that if any time-dependent Hamiltonian be the combination of the generators of any closed 
algebra, an invariant would be a linear combination of the generators of that closed algebra 
with time-dependent coefficients. 

The integrability of the one-dimensional TDHO is not surprising, because this is a lin- 
ear system. We know, one- dimensional time-dependent Hamiltonians usually lead to non- 
integrability, e.g., a simple pendulum whose length varies in time. Except for the adiabatic 
or small oscillation approximations, this was the problem posed by Lorentz at the above 
mentioned Solvay congress and the solutions have the possibility of displaying chaos. 

Now a natural question is if there exist one-dimensional time-dependent nonlinear Hamil- 
tonians which are also integrable? In fact, the singular oscillator with a centrifugal force 
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potential provides an important example whose kinematics is still within the SU(1, 1) alge- 
bra. In the second section, we have discussed such a nonlinear Hamiltonian and we have also 
derived its invariant from our new Lie-algebraic interpretation of the LR- invariant. Using 
this Hamiltonian we have constructed an integrable, area preserving, nonlinear map and 
derived its exact invariant with a knowledge of its fixed points. 

We have determined the classical time-evolution (Perron-Frobenius) operator following a 
recent method elucidated by Rau M. This method was originally used for solving the quan- 
tum Liouville-Bloch equation, but we have applied this in the classical case and have been 
able to show the one-to-one correspondence between this nonlinear map and the linear map 
(derived from the Hamiltonian of the one-dimensional TDHO). In the third section, we have 
studied the quantum dynamics of the nonlinear map and again we have used Rau's method 
to determine the quantum time-evolution operator. This identical approach to classical and 
quantum dynamics helps us to show, rather trivially, the exact one-to-one correspondence 
between the classical and the quantal. Our study also has the dimension of interpreting the 
mathematics of the SU(1, 1) algebra from a dynamical systems viewpoint. For instance, 
we derive powers of non-commuting products of exponentials of the 577(1, 1) algebra. This 
reveals when the powers can degenerate to the identity, something that is immediately clear 
from the dynamics of the underlying system due to the presence of degenerate periodic 
orbits. 



II. THE CLASSICAL KICKED SINGULAR OSCILLATOR 

First, we will describe the general Hamiltonian of which an important special case con- 
stitutes the rest of the paper. For the existence of an invariant for the one-dimensional 
TDHO, the subsidiary condition Eq. (|5|) has to be integrated with some initial conditions, 
that is we should be able to determine a particular solution. The subsidiary condition, given 
in Eq. @, is a nonlinear, time-dependent, equation and its integrability is not immediately 
obvious. If we assume p as a position variable, say q, then Eq. (|5]) is the equation of motion 
corresponding to the Hamiltonian given by, 

1 / o k \ 1 



Hi) = ~lp* + ^)+ -u\t)q\ (6) 
with k = 1. This is a time- dependent Hamiltonian, and has also been studied for long: 



it was studied in part by Lewis and Leach f22|l , Camiz et. al. and Pedrosa et. al. |23 
studied this Hamiltonian quantum mechanically. The new nonlinear force in the system is 
a 'Centrifugal force ' and it appears in many integrable systems, including the celebrated 
Calogero-Sutherland-Moser [p4-26fl many-body Hamiltonian 



The Hamiltonian in Eq. @, is formed by the dynamical variables | (p 2 + and | q 2 . 

These two variables together with pq, also form the SU(1, 1) closed algebra. Following our 
algebraic interpretation of I(t), we can assume the invariant of this nonlinear Hamiltonian 
to be of the form, 

I if) = \a{t) (> + A) + \(3{t)q 2 + 2l(t)pq- (7) 
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Again, by the same substitutions and identical procedures we get 



m = \p 2 [f + + \ [p 2 + j) q 2 - pppq, (8) 

with the same subsidiary condition as given in Eq. (Eh. It is also possible to determine this in- 



variant using time-dependent canonical transformations ||27[| , however, our algebraic method 
accomplishes this more elegantly. In this case, the equation of motion and the subsidiary 
condition are the same nonlinear equation, but the fact that we need only one particular 
solution of this equation to determine the invariant shows the power of the methodology 
and the integrability of the above Hamiltonian in Eq. (|6|). In Fig. [I] we show a special case 
corresponding to to 2 (t) = 1 + cos(\/2t) where the existence of the invariant is reflected in 
the regular structures. 



A. The integrable discrete system 

Now we use the above to construct a nonlinear integrable map. Integrable discrete 
systems, or maps, have been intensely studied for some time now and many interesting 
methods and results have been found. Integrable maps are important for the following 
reasons. Firstly, studying the dynamics of a mapping is more simpler than the dynamics 
of a continuous system since it involves direct iteration. Secondly, if one wants to study 
numerically any integrable continuous system, it is absurd to use non-integrable numerical 
schemes which destroy the basic properties of the system. Therefore, for numerical studies 
of any integrable system, one has to find a discrete integrable version of the system. The 
discrete map we discuss below, or extensions thereof, for instance, may be used in the 
numerical studies of the Calogero-Sutherland-Moser model. Lastly, we can argue that the 
discrete systems are more fundamental ones since they contain continuous ones as special 
limits. 

Studying the stroboscopic map of any time-periodic system is always possible, at least 
numerically. To get an analytical expression a much studied form of the time dependence 
is: oj 2 {t) = uj 2 T ^2 n S(t ~ n T)- The standard map and its quantization as well as an 
experimental realization of the time- dependence through a periodically pulsed laser field 
p8| provides a well known example of this kind. The map corresponding to the TDHO is a 
linear map, which is very simple classically. 

Here we will study the map corresponding to the nonlinear system. The Hamiltonian of 
interest is, 

H(t) = Up 2 + X)+ l^Tq^it-nT), (9) 
and the corresponding Hamilton's equations of motion are, 

P= ^ ~ u 2 T q^25(t - nT), 

^ n 

q = p. (10) 
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The equation of motion of the system is given by, 



q + uj 2 T q^5{t-nT) - A = Q. (11) 



n 



Q 



This equation is the same as that of the subsidiary condition given in Eq. fl5|), except for 
the constant k. Integrating Hamilton's equations of motion from just after the n-th kick to 
just after the (n + l)-th kick and defining new scaled variables 

p^k^T'zp, and q—^k^T^q, 

the phase space map of the system is, 



1 



Qn+l = \\VI + — + Qn + 2 PnQr, 



Pn+i = ^V+i, (12) 

where Q = u T. 

Note that the scaling has removed the k dependence, but that the scaling is singular at 
k = and in fact this limit leads to the removal of the 1/q 2 term, as is discussed further 
below. If we plot this phase space map numerically, it clearly shows its regular behaviour, 
as illustrated in Fig. 0. This rather complicated looking discrete map has a very simple 
behaviour, reflecting the fact that both this and a linear map, derivable from a TDHO, have 
common algebraic antecedents. This linear map is derived from the Hamiltonian in Eq. ([|) 
by simply discarding the 1/q 2 term (the k = limit) and is given by 

q n+ i = p n + q n , 

Pn+l = Pn ~ Q 2 q n +1- (13) 

In fact the property of the linear map that exclusively quasi-periodic or entirely periodic 
behaviour exists for different values of the parameter is also observed for the nonlinear map. 
For < Q 2 < 4 the motion is bounded and stable and for all other values the motion is 
unbounded and unstable, in both the maps. The invariant in the case of the linear map 
describes either an ellipse or a hyperbola, now we determine the invariant of the nonlinear 
map using the method of LR. 

Using Eq. (Rl) the invariant of this map would be of the form, 



P 2 ( Pn 2 + ^2 J + if + ^2 ) <?« 2 ~ ^PPPnQn, (14) 



where p satisfies the same subsidiary condition as given in Eq. (H), but now u 2 (t) = 
u 2 T Yln^it ~ TiT). If we are able to determine any particular solution of the subsidiary 



condition, then we can use that solution to get the invariant of the map Eq. (12). As we 



have already mentioned, the equation of motion of this nonlinear kicked system is the same 



as that of the subsidiary condition for p, as shown in Eq. (|TT|), therefore if we are able to 



get any solution of the nonlinear map Eq. (|12D, that solution should also be the solution 
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P* = — ^174 - (15) 



for p. The most simple solution corresponds to the fixed-point. Therefore the problem of 
determining the invariant has been reduced to determining the fixed-point of the nonlinear 
mapping. The fixed point is at, 

* = 2 

9 fti/2[4(4 - 2 )] 1/4 ' 

[4(4 -W)Y 

This is a particular solution of p and p. Note that unlike the linear map whose fixed point 
at the origin is independent of the system parameters, the fixed points here move with 
the parameter and as the system approaches instability (Q 2 —>■ 4), they approach infinity. 
Using these solutions, we get the invariant of the mapping within an arbitrary multiplicative 
constant as, 

I = p 2 n + \ + n 2 (q 2 n + q nPn ). (16) 

The lack of a time subscript on I indicates its constancy. This invariant is valid even when 
the motion is classically unstable. 



B. The classical evolution operator 



We would like to study the nonlinear map from a point of view that simultaneously 
explains its dynamics as well as sets the stage for quantum mechanical work. Here we 
introduce the classical time-evolution operator for studying the dynamics of the system. For 
a given dynamical variable, say V, the corresponding Liouville operator is denoted by Ly 
and it is defined as Ly = {V, *}, where {*, *} denotes the usual Poisson bracket. 

We know the dynamical equation for any dynamical variable, say /, is given by, 



-L H f, (17) 

where H is the Hamiltonian of the system and Lu is the Liouville operator corresponding 
to H. Let us define A = \{p 2 + £),B = \ q 2 and {A, B} = -pq= -2C, i.e., C = -\pq. 
The triad (A, B, C) form the closed SU(1, 1) algebra. Then we can write, 



L 



H 



L A + u 2 T L B ^5{t-nT) 



(18) 



where La and Lb are the Liouville operators corresponding to A and B respectively. These 
Liouville operators together with Lq, the Liouville operator corresponding to C, forms 
the same 577(1,1) algebra as that of (A,B,C) but under the Lie bracket operation, i.e., 
[L A ,L B ] = -2L C , [L C ,L A ] = L A and [L C ,L B ] = -L B . 

From Eq. fllTD, the dynamical equation for the classical time-evolution of any arbitrary 
function / on the phase space would be, 



dt 



-L H f = - 



L A + u 2 TL B ^5(t-nT) 



(19) 
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Integrating the above equation Eq. fll9|) in between the time t = to t = T, we get the 
classical time-evolution operator from just after zero time to just after the first kick. We 
write f(q(t = T),p(t = T)) = Ff(q(t = 0),p(t = 0)) where 

F = exp (-u 2 T L B ) exp {-L A T) . (20) 

This Perron-Frobenius operator can be regarded as the classical Flouquet operator. To 
understand the dynamics for time n, we have to determine the power F n . The one parameter 
Abelian group of the powers completely specifies the dynamics at all time. However, F is 
itself a product of the exponential of two non-commuting operators which do not even 
commute with their commutator Lq. We can write F in a single exponential form, following 
an earlier work of Truax but that form is so involved that it becomes difficult to extract 



the essentially simple dynamics of the system. 

We apply a recent operator method, referred to earlier to derive the classical time- 
evolution operator at any time. The general procedure of this method is simple to describe. 
If we have any general time-dependent Hamiltonian of the form given by, 

n 

H(t) = ^OityHi, (21) 
i=i 

where (cij(t), i = 1, . . . , n) are a set of linearly independent general complex function of time 
and the dynamical variables (Hi,i = 1, ... , n) are the generators of any n— dimensional 
closed Lie algebra. Corresponding Liouville operators {Ln^i = 1, . . . ,n) would also form 
that same algebra. Then the classical time-evolution operator F(t) can be expressed in the 
product form : 

n 

F(t) = HeMbjit) L H .]. (22) 
i=l 

Therefore, in the classical case we can start with the time-evolution operator of the form : 

F{t) = exp[X(t)L B ] exp\Y(t)L c ] exp[Z(t)L A ], (23) 

where X(t),Y(t) and Z(t) are real functions of time. From the initial condition F(0) = 1, 
we have X(0) = F(0) = Z(0) = 0. Now, substituting this product form of F(t) in Eq. (0), 
and repeatedly applying the Campbell-Baker-Hausdorff (CBH) formula we can cast it into 
a form such that F(t) is pushed to the extreme right in the LHS of the equation Eq. (|TD|). 
This yields a set of first order differential equations for the introduced functions of time : 

X = -X 2 -w 2 T^J(i- nT), 

n 

Y = 2X } (24) 
Z = -e~ Y . 



In the equation Eq. (|2"B|) we can choose the exponential operators in different orders, but 



we found that this leads to sets of differential equations whose solutions may not even exist 
for such kicked systems. 



S 



We now treat the case of delta kicks. Integrating the above equations in between two 
consecutive kicks and defining 

x = TX, y = Y, and z = ^, 

we get a nonlinear mapping, a "coefficient" mapping for the new dimensionless variables 

(x,y,z): 

x n+ i = — — fi 2 , (25a) 

1 + x n 

y n +i =Vn- log [(1 + x n ) 2 ] , (25b) 

z n+l = z n + (x n +l ~ X n ). (25c) 

To be explicit the n-th power of the operator F is 

F n = exp[x n TL B ] exp[y n L c ] exp[z n L A /T]. (26) 

From the initial condition F° = 1, we have Xq = yo = z = 0. The time development is 
now entirely buried in the scalar functions x n , y n , and z n . 

The most important of the recursion equations in Eq. (^) is the first one. We note 
that this equation viewed as a transformation is a special case of the "bilinear" conformal 
transformation in complex analysis or a special case of the projective group PSL(2,M.). We 
solved this nonlinear map by constructing an auxiliary two dimensional linear map. This is 
not entirely surprising as lurking behind the one dimensional nonlinear singular oscillator is 
a two dimensional linear one. This gives a new insight into the often stated close relationship 
between the harmonic oscillator and the singular oscillator. 

Before we solve these equations explicitly we point out the following interesting fact. 
Consider the Hamiltonian 

H(q,p,t) = qp 2 + fi 2 g^5(t-n), (27) 

n 

written in terms of dimensionless canonical variables. The triad (qp 2 , q, qp) form the alge- 
bra 5*0(2, 1) under the Poisson bracket operation and this algebra is locally isomorphic to 



SU(1, 1) p0|j31| . Then the resulting map for q n and p n is 



Pn+i = - n 2 , (28a) 

1 +Pn 

q n+1 = q n (l + Pn) 2 ■ (28b) 

Thus if we identify x n = p n and y n = — \og(q n ) we get the first two of the recursion relations 
in Eq. (^5|). We note that the SO(2, 1) algebra underlies the Morse oscillator p2| |. Therefore 



the coefficient map also has an algebraic basis and in fact one that is practically identical 
to the original one. 

We return now to solve Eq. (^Sj). First we define s n = 1 + x n , and in terms of s n , Eq. 
( |25a|) becomes 

1 . . 

s n +i = V , (29) 
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where rj = 2 — Q 2 , and from the initial condition xq = 0, we have so = 1. Construct the 
auxiliary linear map: 



( a-n+i \ = ( V -1 \ ( a r , 
V b n+1 J \ 1 J V K 



(30) 



with initial conditions ao = bo — 1. We identify s n = a n /b n . 

To get the general form of s n , we have to diagonalize the matrix M = ((?/, —1), (1,0)). 
The eigenvalues of M are, 

A ± = i(r7±vV-4). (31) 

Whether A± is real or complex is dependent on rj. Therefore just as for the linear map we 
have to study separately three different regions for the parameter rj, these are : 

Case 1: rf < 4, i.e., < Q 2 < 4. 

In this case the eigenvalues of M are complex-conjugate of each other, i.e., we can write 

(1-Vl 2 

A± = exp(±«cr), where a = cos 1 (r]/2) = cos 1 lf 



The dynamics is stable and bounded. After diagonalizing the matrix M, we determine a n 
and b n to finally get: 



Xn = cos a — sin a tan 



2 

In terms of x we can obtain the solution for y and z as : 

n-l 



( 2n - 1)^-1. (32) 



y n = -2^1og|l , (33) 

k=0 

= (34) 

When a = 2iTm/N, where m and N are coprime integers, one sees, after some algebra, 
that 

X n = Xn+N, Dn = Dn+N-, Z n = ^ n +V- 

Therefore for these particular values of o (or of the corresponding value of Q), the above 
three-dimensional map for the coefficients is exactly periodic. To be explicit at these values 
of the parameter 

F N = 1, (35) 

the Perron-Frobenius operator becomes unity at time N. The quantum equivalence to be 
discussed below will be complete and the Flouquet unitary operator will become unity at 
the same time. Fig. ^ shows the coefficients for a periodic case. Since the coefficient z is 
proportional to x, it would follow the same behaviour as x. This map in Eq. (O) would also 
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be periodic with the same period, and corresponds to the phase space shown in Fig. @(a). 
For other values of Q, this map is quasi-periodic. Equivalently, the phase space map Eq. 
(|T2|) also displays the above behaviour for corresponding values of Q. This shows that our 
operator approach for studying the classical map is in one-to-one correspondence with the 
phase space dynamics. Fig. |] shows the coefficients for a quasi-periodic case, corresponding 
to the phase space in Fig. |2](b). Again the behaviour of z is same as x. 

Case 2: rf > 4, i.e., Q 2 > 4 or Q 2 < 0. 

We can divide this case into two parts. They are : 

(a) 7] > 2, i.e., Q 2 < : 

In this part the eigenvalues of M are real and positive. Therefore we can take 

\± = exp(±cr), where a = cosh 1 (r]/2) = cosh 1 ( — - — j . 

In this part the dynamics is unstable and unbounded (hyperbolic). Again following proce- 
dures as outlined above, we get 



(2n-l)|] - 1. (36) 



x n = cosh a + sinh a tanh 

In this case the solution for y and z would also be the same as that given in Eq. (|33]). How- 
ever, the basic properties of x, y and z would change due to the unstable and unbounded 
dynamics, and this is evident from the above equations. For large n, both x and z asymp- 
totically reach a constant value that depends on the magnitude of Q. 

However, y would increase linearly with n, when n is large, 
(b) 7] < -2, i.e., tt 2 > 4. 

Here both the eigenvalues of M are real and negative. Thus we can take 

A± = — exp(±cr) where a = cosh~ 1 (|/7|/2) = cosh -1 I j . 

In this part the dynamics is still unstable and unbounded; it corresponds, in the linear map, 
to hyperbolic fixed point with reflections. On following the above procedures, we get 



x„ = — cosh a — sinh a coth 



(2n-l)|l - 1. (37) 



The solutions for y and z still remain the same, and their behaviour for large values of n are 
qualitatively same as that of the previous part. 

Case 3: r] 2 = 4, i.e., H 2 = or fi 2 = 4. 

These cases correspond to the marginal ones separating the stable and unstable motions. 
We can divide this case also into two parts. They are : 
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(a) r] = 2, i.e., ft 2 = 0. 

Here the eigenvalues of M are equal, A± = 1 . The kick is not operating on the system. This 
implies that in the expression of the time-evolution operator F, x = y = 0. Therefore F 
contains only one exponential and hence the coefficient z would increase linearly with time. 
This can also be seen as a limiting case; as from Eq. ([25 aj) , (x n +i — x n ) — > —Q 2 and the Eq. 
gives z n = —n. 



(b) 7] = -2, i.e., VI 2 = 4. 

Again, the eigenvalues of M are equal, but now A± = —1. We can get the solution for the 
coefficients quite easily and these are given by, 

2n + 1 . . 

x n = 7 ~ 1) (38a) 

2n — 1 

y n = -2 1og|2n-l|, (38b) 

z n = (38c) 

For large n, the coefficients x and z asymptotically reach constant values (—2 and —1/2), 
while the magnitude of y increases logarithmically. Thus this marginal case straddling the 
stable and the reflective hyperbolic cases would have power law behaviors in time for phase 
space variables. 



III. QUANTUM DYNAMICS OF THE KICKED SINGULAR OSCILLATOR 

Quantum mechanical studies of the time-dependent singular oscillator have been carried 
out for some time now, for instance in |23] where complete analytical solutions were given. 
We exploit our algebraic method to the special case of the kicked oscillator to lay bare 
properties such as exact periodicity and quasi-periodicity. In fact, except for a change of 
terminology, the mathematics is already complete in the previous section. 

During our study of the classical nonlinear map, we had introduced the classical time- 
evolution operator to show the one-to-one correspondence between the nonlinear and the 
linear map. While classically this is not the usual approach to dynamics, in the case of quan- 
tum dynamics, the most natural and popular way is to study the quantum time-evolution 
operator U (£). Thus our classical approach generalizes most easily to the quantum. Previous 
work |33[ that points out exact quantum-classical correspondence in the case of the SU(1, 1) 
algebra for the coefficients of the invariant is easily understood in our approach. 

We define the operators as in the classical case: 

A = \ if + ' ^ = 2^' = ~\ h{S>( l + 9P) = 

where (A, B, C) form the closed SU(1, 1) algebra which is given by [A, B\ = 2HC, [C, A] = 
hA and [C, B] = —KB. Among these operators C is anti-hermitian and the other two are 
hermitian. In terms of these operators our Hamiltonian would be, 

H(t) = A + u 2 TB ^8{t-nT), (39) 

n 
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The time-evolution operator U (t) satisfies the equation, 



at 



(40) 



where £7(0) = 1 and U(t) is unitary. The quantum Flouquet operator U(t = T), which we 
denote simply as U, is given by 



U = exp ( —-u TBj exp I - J.T A 



(41) 



For the complete dynamics, as usual, we have to determine the powers U n . Again we apply 
Rau's IJ method for the derivation of the time-evolution operator at any arbitrary time. 
We start with the time-evolution operator of the form, 



U(t) = exp 



- n X (t) B 



exp 



- h Y{t)C 



exp 



-Z (t )A 



(42) 



where X(t),Y(t) and Z(t) are real functions of time, so that U(t) remains unitary. From 
the initial condition U(0) = 1, we have X(0) = Y(0) = Z(0) = 0. Substituting the above 
U(t) in Eq. fl40|), following identical procedures as given for the classical dynamics, we get 



X = -X 2 -iu 2 T^25(t- nT), 

n 

Y = 2X, 



(43) 



Z = -e 



-Y 



These set of equations for the quantum coefficients are identical to the equations for the 
classical coefficients as given in Eq. (p4[). Therefore we would get an identical map for the 
quantum coefficients as for the classical coefficients Eq. fl2"5]). Since these two maps are 
identical, their stability properties are identical, i.e., quantum dynamics exactly follows its 
classical counterpart. Thus when the classical Perron- Frobenius operator becomes identity 



as in Eq. (|35|) the Flouquet operator also becomes identity. The various cases discussed 
classically, including the marginal and the reflective hyperbolic cases have exact quantum 
counterparts. Thus as far as time evolution is concerned the quantal problem is already 
solved. Issues of eigenstates and spectrum can be tackled similarly, as for example done in 
3~4| for the case of the kicked harmonic oscillator; however we do not pursue this further 
here. 



IV. SUMMARY 

We have studied, both classically and quantum mechanically, a time- dependent one- 
dimensional nonlinear integrable system, namely the kicked harmonic oscillator with a sin- 
gular potential. By applying our Lie-algebraic interpretation of the LR-invariant, we have 
determined the exact invariant of that nonlinear system. One can apply this method to de- 
termine the invariant of any time-dependent Hamiltonian, formed by the generators of any 
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closed algebra. We have also found, as far as we know, a new integrable nonlinear mapping. 
Here we derived that map from the Hamiltonian formed of the generators of the SU(1, 1) 
algebra. We expect that we can determine integrable two-dimensional mappings from the 
Hamiltonians formed by the generators of other closed algebras. 

We constructed the classical time-evolution operator, or the Perron-Frobenius operator 
for the nonlinear integrable system, taking advantage of a recent method of solving the 
time-dependent Schrodinger equation. This brings into relief the extreme quantum-classical 
correspondence in these systems which have algebraic structures underlying them. 

We may speculate whether for time-dependent systems with more than one degree of 
freedom, which are constructed from elements of some Lie algebra, the nonlinear equations 
(whether difference or differential) that are the analogue of the coefficient mapping in Eq. 
(E5|) are integrable. Further, if there is a connection between the integrability of these 
equations and that of the original system itself. Our study that began as an attempt to 
study a nonlinear time- dependent integrable system led us to one that is in many respects 
related to the harmonic oscillator with time-dependent frequency. We may then ask if we 
can go beyond this limitation to a more "genuine" form of nonlinearity. 
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FIGURES 



FIG. 1. Stroboscopic picture of the time-dependent harmonic oscillator with a singular pertur- 
bation; the integrable behaviour is evident. All quantities shown are dimensionless 

FIG. 2. (a) Period-9 orbit of the nonlinear map for the case n 2 = 2 [i — cos (^p)j ~ 0.4679, 
(b) quasi-periodic orbit for the case O 2 = 2.43. All quantities shown are dimensionless 

FIG. 3. The coefficients of the time-evolution operator showing the period-9 behaviour for the 
same O 2 as in Fig. §(a). Shown are the x (a) and y (b) coefficients. All quantities shown are 
dimensionless 

FIG. 4. The coefficients of the time-evolution operator showing the quasi-periodic behaviour 
for the same Q 2 as in Fig. |2|(b). Shown are the x (a) and y (b) coefficients. All quantities shown 
are dimensionless 
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